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Weak solution for compressible fluid models of Korteweg 

type 

Boris Haspot * 



Abstract 



This work is devoted to proving existence of global weak solutions for a general isothermal 
model of capillary fluids derived by J.E Dunn and J.Serrin (1985) [5J, which can be used 
as a phase transition model. 

We distinguish two cases when the dimension N = 2 and N = 1, in the first case we need 



that 1 G L°°, when N = 1 we get a weak solution in finite time in the energy space. 

1 Introduction 

1.1 Derivation of Korteweg model 



We are concerned with compressible fluids endowed with internal capillarity. The model 
we consider originates from the XlXth century work by van der Waals and Korteweg 
|1U| and was actually derived in its modern form in the 1980s using the second gradient 
theory, see for instance (9j [E] . 
Korteweg-type models are based on an extended version of nonequilibrium thermodynam- 
ics, which assumes that the energy of the fluid not only depends on standard variables 
but also on the gradient of the density. 

The model derives from a Cahn-Hilliard like free energy (see the pioneering work by 
J.E.Dunn and J.Serrin in [5] and also in [HOE]), the conservation of mass reads: 



CO 

o 

OO 

O 



and: 



d 

—p + div(pu) = 0, 

(1.1) 

d 

— {pu) + div (pu (g) u) — pAu — (A + p)Vdivu + VP(p) = npVAp, 
*. ot 



p > and 2p + A > 

where p represents the density, u G M. N the velocity, P is a general pressure, p, A are the 
coefficients of viscosity and k the coefficient of capillarity. The term npVAp corresponds 
to the capillarity term and allows to describe the variation of density in the interfaces 
between two phases. 
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One can now rewrite pVA/9 in the following form to understand the difficulty of the non 
linear terms in distribution sense: 



KVAp = divK, 



with Kij = -^(Ap 2 - \Vp\ 2 )5ij - Kdipdjp. 



X2) 



One can recall the classical energy inequality for the Korteweg system. Let p > be a 
constant reference density, and let II be defined by: 



so that P(s) = sU'(s) - n(s) , Il'(p) = and: 

d t U(p) + div(ull(p)) + P(p)div(«) = in £>'((0, T) x R N ). 

Notice that IT is convex as far as P is non decreasing (since P (s) = sll"(s)), which 
is the case for 7-type pressure laws or for Van der Waals law above the critical tempera- 
ture. 

Multiplying the equation of momentum conservation in the system (jl.ip by pu and inte- 
grating by parts over Mr, we obtain the following estimate: 

/ [\p\u\ 2 + (U(p) - U(p)) + ^\Vp\ 2 )(t)dx + 2 f I {2pD(u) : D(u) 

JR N 1 1 JO JR N 

+ (A + p)\divu\ 2 )dx < [ + (n(po) - n{p)) + ^\v P0 \ 2 )dx. (L3) 

Jrjv Zp I 

We will note in the sequel: 

£(t)=[ C-p\n\ 2 + (U(p)-U(p)) + ^\Vp\ 2 )(t)dx, (1.4) 
J r n v 2 2 

It follows that assuming that the initial total energy is finite: 
r l 1 2 

eo = / (M- + (n(p ) - n(p)) + ^\Vp \ 2 )dx < +00 , 
Jrn 2p 2 

then we have the a priori bounds: 

n(p)-n(p), and p\u\ 2 e L 1 (0,oo,L 1 (R N )), 

Vp€ L oo (0,oo,L 2 (R N )) N , and Vu G L 2 (0,oo,M. n ) n2 . 

In the sequel, we aim at solving the problem of global existence of weak solution for the 
system (jl.ip . Assuming that we dispose from a smooth approximates sequel (p n ,u n ) n€ fq 
of system one can remark easily that it is difficult to pass to the limit in the 
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quadratic term V p n <g> V p n which belongs to L 00 ^ 1 ). According to the classical the- 
orems on weak topology, V p n <S> Vp n converges up to extraction to a measure u, the 
difficulty is to prove that v = V p V p where p is a limit of the sequence (p n )neN in 
appropriate space. 

Another difficulty in compressible fluid mechanics is to deal with the vacuum and we will 
see that this problem does appear in the model of Korteweg, when estimating V/>. As a 
matter of fact, the existence of global solution in time for the model of Korteweg is still 
an open problem. 

The first ones to have studied the problem, are R. Danchin and B. Desjardins in [1]. 
They showed that if we take initial data close to a stable equilibrium in the energy space 
and assume that we control the vacuum and the norm L°° of the density p, then we get 
some weak solution globally in time. Controlling the vacuum amounts to bounding = in 
L°°. 

Recently D. Bresch, B. Desjardins and C-K. Lin in [2] got some global weak solutions 
for the isotherm Korteweg model with some specific viscosity coefficients. In effect, they 
assume that p(p) = Cp with C > and X(p) = 0. In choosing these specific coefficients 
they can get a gain of derivative on the density p and obtain an estimate for p in L 2 (H 2 ). 
It is easy now with this type of estimate on p to pass to the limit in the term of capil- 
larity. However a new difficulty appears with the loss of information on the velocity u 
and it becomes difficult to pass to the limit in the term pu (8> u. Hence the solutions of 
D. Bresch, B. Desjardins and C-K. Lin require some specific test functions which depend 
on the density p. Indeed the loss of information is on the vacuum. 

Our result is in the same spirit as the one by Danchin and Desjardins: we want to 
improve the energy inequality which allows a gain of derivative on the density p. We 
show that we don't need to control p in L°° norm to get global weak solution. However 
the control of the vacuum seems necessary. 

In section [2] we recall some definitions on the Orlicz space and some energy inequality 
on the system in these spaces. In section [3] we show a theorem on the global existence of 
weak solutions in dimension two under some conditions which amount to controlling the 
vacuum. In the last section we investigate the case of the dimension one, and we get a 
theorem of local existence of solution in the energy space and a result of global existence 
with small initial data. 

2 Classical a priori estimates and Orlicz spaces 
2.1 Classical a priori estimates 

We first want to explain how it is possible to obtain natural a priori bounds which 
correspond to energy when the density is close to a constant state. 

We first rewrite the mass equation in using renormalized solutions, and the momentum 
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equation. We get the following formal identities: 



7 — 1 Ot v ' ' ' L 7 

a M 2 id 2 ( 2 - 5 ) 

p— — - — h pu • V— pAu • u — £Vdivii • u + au • Vp 7 = Kp-u • VAp , 



where we note: £ = p + A. 

Therefore we find in summing the two equalities of (12, 5ft : 



I + ;^-[(P 7 + (7- l)/^ 7 " 7P 7 ~V)] + div^-^p 7 - p^p) + P^)] 

— pAu • u — £Vdivu • u = KpVAp • u. 

(2.6) 

We may then integrate in space the equality (|2.6[) and we get: 



u\ dx 



I {p^ + ^{p^ + {l-l)p 1 -lfP- 1 P ) + K\Vp\ 2 )(t,x)dx + pf ds f \D 
Jr n 27-1 J Q J R N 

+ tfdsf \&ivu\ 2 dx < [ (p ^ + -%( P 2 + (7 - 1)P 7 - J^Po) 
Jo Jr n Jr n 2 7 - 1 



+ k\ Vpo| 2 ) (x) dx. 

(2.7) 



Notation 1 In the sequel we will note: 

3j(p) = P 1 + (7 - 1)P 7 - 7P 7 ~V 



\u\ 2 a 



£ 7 (t)= / (p^ + — ( p y + ( 7 -i)fF- 1 p / -i p ) + K \v P \ 2 )(t,x)dx 
Jrk 2 7-1 

^ = / {P^f- + ^7W + (7- 1)^ " 7P 7 - 1 Po)^ 
Jr^ 2 7 — 1 



We now want to estimate this quantity j-y(p) and in this goal we recall some properties 
of Orlicz spaces. 

2.2 Orlicz spaces 

We begin by describing the Orlicz space in which we will work: 

L%R N ) = {/ e Lj oc (R N )/n { \f\<s} e LF(R N ), fl {lfm e L 9 (R N )} 
where 5 is fixed, 5 > 0. 

First of all, it is not difficult to check that Lp does not depend on the choice of 5 > since 
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is bounded from above and from below on any interval [Si, 62] with < Si < S2 < +00. 
In particular we deduce that we have: 

f G Ll(WL N ) if / € L|(M JV ') and p, q > e. 

s 

Obviously we get meas{|/| > S} < +00 if / G Lp(R N ) and thus we have the embedding: 

L q p (R N ) C L q p \(R N ) if 1 < qi < q < +00, 1 < p < Pl < +00. 

Next, we choose ^ a convex function on [0, +00) which is equal (or equivalent) to x p for 
x small and to x q for x large, then we can define the space Lp(R N ) as follows: 

Definition 2.1 VFe define then the Orlicz space Lp(M. N ) as follows: 

L q (R N ) = {/ G LL(M iV )/^(/) G L 1 ^)}- 

We can check that Lp(IR ) is a linear vector space. Now we endow L q p (R N ) with a norm 
so that L^R-^) is a separable Banach space: 

ll/llL^)=mf{t>0/ *({)<1}. 

We recall now some useful properties of Orlicz spaces. 
Proposition 2.1 The following properties hold: 

1. Dual space: If p > 1 and q > 1 then (L q p (R N ))' = L q ,(R N ) where q = p = 

p 

P-i- 

2. L q = L p + L q ifl<q<p< +00. 

3. Composition: Let F be a continuous function on R such that F(Q) = 0, F is 
differentiate at and F(t)\t\~ e —> a / at t —> +00. Then if q> 9, 

F(f)e4(R N ) if feLl(R N ). 

Now we can recall a theorem on the Orlicz space concerning the inequality of energy 
Proposition 2.2 The function j 1 {p) is in L 1 (R Ar ) if and only if p — p G l? 2 . 
Proof: 

On the set {\p — p\ < S}, p is bounded from above, since 7 > 1 we thus deduce that j 7 (/c) 
is equivalent to \p — p\ 2 on the set {\p — p\ < S}. 

Next on the set {\p — p\ > S}, we observe that for some v G (0, 1) and C G (1, +00), we 
have: 

^P-P^ <h( P )<C\p-p\\ 

□ 
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Link with our energy estimate 

We recall the definition of the fractional derivative operator A s : 
Definition 2.2 We define the operator A s as follows: 

We recall now some useful results, we start with a proposition coming from the theorem 
of interpolation by Riesz-Thorin. 

Proposition 2.3 The Fourier transform is continuous from LP in L q with p G [1,2], 
q G [2, +oo] and: 

1 1 

- + - = 1. 

P Q 

We recall the definition of homogeneous Sobolev space. 
Definition 2.3 Let s£R. f is in the homogeneous space H s if: 

\tff€L 2 (R N ). 

Proposition 2.4 Let f G H s with s > and f G L p + L 2 with 1 < p < 2. Then f G L 2 . 
Proof: 

Indeed we have as / G H s : 



[ \e s \f\ 2 da < +oo 

Jr n 



so j| >1} 6 L 2 ^). And as / = /i + / a with /i G LP(M 7V ) and / 2 G L 2 . In using the 
Riesz-Thorin theorem, we know that /i G L 8 ^) with 1 + 1 = 1. 



As g > 2 we then have /l{i/j<n ^ (M ). This concludes the proof. □ 

According to the above theorem and our energy estimate we get that for all T G R, 
p-peL°°(p,T;L2 



Remark 1 We have then in using previous properties on Orlicz spaces and (2.7): 
• i/7 > 2 i/iera Lj(R^) <-► L 2 (M Ar ) and so p — p £ L°°(H 1 (R N )). 



• 7 < 2 i/ien following the proposition 2.4 and the fact that LS^ = L 7 + L 2 we <?e£ 
p-p€ L 00 (^ 1 (M JV )). 

We can now explain what we mean by weak solution of problem (jl.ip in dimension 
N = 1,2. 

Definition 2.4 Let the couple (pq,uq) satisfy; 

1. po G Ll(R N ), V Po G L 2 (R N ) and ± G L°°(]R iV ). 



2. p Q \u \ 2 G L\R N ) 

3. PqUq = whenever x G {po = 0}, 
We have the following definition: 

• A couple (p, u) is called a weak solution of problem on I x R N with I an 
interval ofR if: 

- p G L°°(Ll(R N )), Vp G L°°(L 2 (R N )), <pp G L 2 (iJ 1+a (R Ar )) Va e]0,l[, and 
V99 G Cg°(R Ar ). 

- i G L°°((0,+oo) x R^). 

- Vn G L 2 (L 2 (R^)), p\u\ 2 G L°°(L l ). 

- Mass equation holds in T> (I x M N ). 

- Momentum equation holds in T> (I x M. N ) N . 

- Hm t _ + J RN p(t)<p = J RN p (f, W(f G V(R N ), 

- lim t _ >0+ f RN pu(t) ■ 4> = f RN {pu)o • <p, y<t> G V(R N ) N . 

• The quantity £q is finite and inequality \2. 7\ ) holds a.e in I. 

3 Existence of global weak solutions for N = 2 
3.1 Gain of derivability in the case N = 2 

In the following theorem we are interested by getting a gain of derivative on the density 
p. This will enable us to treat in distribution sense the quadratic term V/9 (g> Vp. 

Theorem 3.1 Let N = 2 and (p,u) be a smooth approximate solution of the system 
(jj.jp such that - G L°°((0, T) x R N ). Then there exists a constant 77 > depending only 
on the constant intervening in the Sobolev embedding such that if: 

I|Vpo||l 2 (r 2 ) + II\/PoKIHl 2 (i; 2 ) + IIj'tCpo)!^ 1 < V 
then we get for all ip G C^{R N ): 

||</>p 2 || i2(H i +4) < M with < s < 2, 

where M depends only on the initial conditions data, on T, on ip, on s and on \\-\\l<x>. 

Remark 2 In fact instead of supposing that -, we can just assume that Vlogp G 
L°°(L 2 ). This will imply that p will be a weight of Muckenhoupt. 

Remark 3 In the sequel the notation of space follows those by Runst, Sickel in 
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Proof of theorem I3.lt 



Our goal is to get a gain of derivative on the density in using energy inequalities and in 
taking advantage of the term of capillarity. We need to localize the argument to control 
the low frequencies. Let ip E C^°(M. N ) and we have then: 

d t div((ppu)+di } j(ippuiUj) - (A + 2p)A(ipdivu) + A(<pP(p)) 

= (A(^ 2 ) - ^|Vp| 2 ) - Kdl^di^p)) + (3 ' 8) 

with: 

d 

Rtp = ~q^{p u ' ^V) + \^i,jf)P u i u j + 2di,(pdj(puiUj) — (2p + A)A<pdivu 

- 2(2/i + X)Vtp ■ Vdivu + Ai^ap 7 + 2aVip ■ V(p 7 ) - ^Atp(Ap 2 - \Vp\ 2 ) 

- KVip ■ V{Ap 2 - \Vp\ 2 ) + K(d£ d (p)dipdjp + 2Kdi(pdj(dipdjp) - kA(V^ • Vp 2 ) 

We can apply to the momentum equation the operator A _1 (A) _1 div in order to make 
appear a term in Ap 2 coming from the capillarity. Then we obtain: 



^A( W 2 ) + ^A-%|Vp| 2 ) + Kk^RiRiiipdipdjp) = -A-^div^) 
2 2 dt (3.9) 

+ A^RiRjfapuiUj) - (2p + A)A- 1 ( ( pdivu) + A _1 (^P(p)) + A^fA) -1 ^, 

where Ri denotes the classical Riesz operator. We multiply now the previous equality by 
A 1+S (cpp 2 ) and we integrate on space and in time: 

£ f I \k 1+ Hvp 2 )?dxdt + ^- [ T I {\^p\ 2 )K s ^p 2 )dxdt 
1 Jo Jr n 1 Jo Jr n 



[ V R i R j (^d lP d j p)A s (ipp 2 )dxdt = [ A- 3 div(<ppu)A 1+s (<pp 2 )(T)dx 
Jr n Tf ' ' Jr n 

A~ 3 div(<pp u )A 1+s (ippl)dx - [ [ A~' 3 div (tp p u) A 1+ s -^(ipp 2 )dxdt 

Jo Jrn ' dt 

+ {2p + \) I ( (pdivuA s (<pp 2 )dxdt- [ [ RiR j ((ppu i u j )A s (<pp 2 )dxdt 
Jo Jr n Jo Jr n ~i 

+ f [ <pP(p)A s ((pp 2 )dxdt+ [ [ {A)- 1 R 9 A s {ipp 2 )dxdt. 
Jo Jr n Jo Jr n 

(3.10) 

Now we want to control the term f Q f RN |A 1+ 2 ((pp 2 )\ 2 . Before getting in the heart of the 
proof we want to rewrite the inequality (|3.10p in particular the term: 

A- 3 div( WU )A 1+s ^(^ 2 ). 



Jr n 
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In this goal we recall the renormalized equation for ipp 2 : 

d 

— (ipp 2 ) + div((pp 2 u) = -(pp 2 divu + r v , (3-H) 



with r<p = —V(p • p 2 u. 

So in using the renormalized equation (13. lip we have: 



[ [ A 3 div (p p u) A 1+s (p p 2 ) dxdt = - [ [ A 2 : div (p p u) A s (pp 2 d\vu)dxdt 
Jo Jrn at J J r n 



T r rT 



A 2 d\\ (p p u) A s div (p p 2 u) dxdt + A 3 di\(<ppu)r ip . 

lo Jr n Jo Jr n 

(3.12) 

In combining (|3.1U|) and (|3.12p we get: 

£ / T / |A 1+ i( W 2 )| 2 ^ + ~ / T / (|^Vp| 2 )A s ( W 2 )^di 
z io Jr n 1 Jo Jr n 



K 



r 



o Jr" y 



J2RiR j (ipd i pd j p)A s (ipp 2 )dxdt= / A- 3 div( W n)A 1+s ( w 2 )(T)(ix 



A 3 div(vvo«o)A 1+s (^/9o)cfo- / / A 2 div(^p«)A s ((/?,9 2 divu)<ix(ft 

/ A _2 div((/?pu)A s div(99/) 2 n)(ix(it + (2/i + A) / / pdivuA s (pp 2 )dxdt 
o Jr n ' Jo Jr n 

r rT 



+ [ [ R i R j (ppu i u j )A s (pp 2 )dxdt + [ [ p>P(p)A s (pp 2 )dxdt 
Jo Jr n ~1 ' Jo Jr n 

+ I I (A)- 1 R v A 1+ i(i P p 2 )dxdt+ [ [ A- 3 div(ppu)A l+s r v dxdt. 
Jo Jr n Jo Jr n 

(3.13) 

In order to control J RN | A 1+ 2 (pp 2 )\ 2 , it suffices to bound all the other terms of (|3.13|) . 
Next we will have a control of A 1+ 2 (ipp 2 ) and so a gain of I derivative on the density. 
We start with the the most complicated term which requires a control of - in L°°. 



l)/o T / K -MVp| 2 )A%p 2 ): 

By induction we have V(pp 2 ) G L\(H^-). So by Sobolev embedding we get V(pp 2 ) E 
L 2 (L P ) with = 2 — f (we remark that the case s = 2 is critical for Sobolev embedding) 

Now we have pV ' p = V ^ - — \p\?p and we recall that by hypothesis ~ E L°°, so we 
have ipVp E L 2 (L P ) because pVp E L°°(L r ) for all 1 < r < +00 as p - p E L°°{H l ). 
We now consider A s (pp 2 ). We have by induction A s (pp 2 ) E L 2 (H l ~z) and A s (pp 2 ) G 
L 2 (L 2 ) because I/?/? 2 G L 2 (L 2 ) which enables us to control the low frequencies of A s (pp 2 ). 
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We have then A%p 2 ) G L 2 {H X -*) 
Finally by Holder inequality we get tp\ Vp| 2 A s (ipp 2 ) G L^(L 1 (R N )) because \ + \ + \ = 
\ + \ — | + | = 1 and we have: 

f {v\Vp\ 2 W{w 2 )dxdt < \\V P \\l ¥ (l^%^ 2 )\\l^)\\^p\\\l^l P ), 
Jm (3.14) 

We proceed similarly for the term: 

f [ ^ RiRj (^dipdj P )A s {ipp 2 )dxdt 
Jo 

because we have in following the same lines tpdipdjp G L 2 {L q ) with | = 1 — | and we 
have the Riesz operator which is continuous from L p in L p for 1 < p < +oo. 

We next study the term J RN A~ 3 div((pp u)A 1+s ((pp 2 )(t)dxdt. 
2) f RN A- 3 dw(vpu)A 1+s {ipp 2 )dx: 

We rewrite the term J RN A _3 div(p u)A 1+s (<pp 2 ) on the form: 

f A- 3 div(ppu)A 1+s (vp 2 )dx = [ hT x &v{tp P u)hT 1+8 {ipp i )dx 
Jr n Jr n 

= V / Ri&pu^A-^&p^dx. 

l<i<N J ® N 

As J e Lf{L°°) then we have u G Lf{L 2 ). We recall that p - p G L°°(i? 1 ) then 
<pp G L°°{L P ) for all 1 < p < +oo. We deduce that (pp-u belongs to L°°(L 2 ~ /3 n L 1 ) for 
/? > 0. So we have Ri((ppiii) G L^{L r ) for all 1 < r < 2 by continuity of the operator i?j 
from L p to L p when 1 < p < +oo. 

Case 1 < s < 2: 
Next we have: 

V(<pp 2 ) = 2<ppVp + p 2 Vip 

then we get V(<pp 2 ) G L°°(L 2 ~ 13 ), in using the fact that p — p G L°°(H l ) and Sobolev 
embedding with Holder inequalities. We get then that ipp 2 belongs to L 0O (W 2 1 _ /3 ) = 



We have then A s ~ 1 ((pp 2 ) G L°°(H%_p). By Sobolev embedding we get A 3 " 1 ^ 2 ) G 

L°°(L P ) with | = — 2 +/ £ = 2^/3 + f ~~ ft sma U enough to avoid critical 

embedding. 

Finally we get i^(<pp-u)A- 1+s (<pp 2 ) G L\,{L 1 {E. N ). Indeed we have J + ^ = g^+f-l < 
1 in choosing /3 small enough and ^ + = — 1 + jzp + | > 1 in choosing (3 small 
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enough if necessary, we conclude by interpolation. 
We have finally: 

A- 3 div(tppu)A 1+s (ipp 2 )dx\ < M 



with Mq depending only on the initial data. 
Case < s < 1: 

In this case we conclude by interpolation with the previous case. We now want to study 
the other terms coming from the renormalised equation (|3.11|) . 

3) J Q T J RN A~ 2 div(ippu)A s ((pp 2 divu)dxdt and J RN A~ 2 div(ippu)A s (div ((pp 2 u))dxdt: 
We start with: 

[ A^ 1 div((ppu)A s ~ 1 (div(ipp 2 u)) = [ [ div(ippu)A s ~ 2 (div (ipp 2 u)). 
o Jr n Jo Jr n 

Case 1 < s < 2: 

We have: 

d\v{ippu) = u ■ V(ipp) + (ppdivu. 

By Holder inequalities and Sobolev embedding we get that div(ippu) belongs to L^(L 2 ~^) 
for all p e] 0,1]. 

Next we rewrite div(ipp 2 u) on the form: 

div((pp 2 u) = u ■ V((pp 2 ) + ipp 2 divu. 

As previously dw((pp 2 u) is in L\{I?~^) for all (3 G]0, 1]. Now by Sobolev embedding we 
have A s_2 div(93 / o 2 u) G L^L^) with = = i^rp — 2 -j^ with (3 small enough to avoid critical 
Sobolev embedding. 

We conclude that div(ippu)A s ~ 2 (div((pp 2 u)) is in L^L 1 ) because 2=/J + § = 2^/3 — = 
2^rp — 1 + | < 1 with (3 small enough if necessary and 1 + ^ > 1, so we obtain the result 
by interpolation. Finally we have: 

T r 

/ A~ 2 div(^pit)A s (div(^p 2 «))dxdi| < M 
o Jr n 

with Mq depending only on the initial data. 
Case < s < 1: 

We have the result by interpolation with the previous case. 

Next we proceed similarly for: 

,T r- 

I A 1 div((ppu)A s 1 (tpp 2 d\vu)dxdt. 
o Jr n 
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4) Last terms 

We now want to concentrate us on the following term: 

rT r 

/ RiRj((ppuiUj)A s ((pp 2 )dxdt 

We know that u G L°°(L 2 ) and Du G L 2 (L 2 ) then « G L^il 1 ) and by Holder inequalities 
and Sobolev embedding we can show that ippuiUj G L\{L 2 ) and so RiRj(tppuiUj) G 
4(/- 2 )- 

We have seen that A s ((pp 2 ) G L 2 (i7 1 2 ) then we have as 1 — | > 0: 



|A s (w 2 )ll^(L 2 )<Af + ||^ 2||; 



with < (3 < 1. 
We have then: 



r 

iR^ 



Y,RiRi^pUiU j )K s {ipp 2 )dxdt\ < M + ||A 1+ 5( W 2 )||^ 2(L2) 



with < j3 < 1 and Mo depending only on the initial data. 
We are interested by the term: 

/ tpdivu A s (ipp 2 )dxdt 
Jr n 

We have then diva G L 2 (L 2 ) and we have shown that A s ((pp 2 ) G L°°(L 2 ) so we conclude 
in the same way than the previous term. 
We finally conclude with the term: 

I (pP(p)A s ((pp 2 )dxdt. 
Jr n 

Similarly we have A s (ipp 2 ) G L\-,{L 2 ) and <pP(p) G L\-,{L 2 ) because ipp is in L°°(ff 1 ), and 
we conclude by Sobolev embedding. 

To finish we have to control the term with R^ and r v that we leave to the reader. These 
terms are easy because they are more regular than the preceding terms. We finally get 
in using all the previous inequalities: 

l|A 1+ f (w> 2 )lH 2(i2) < Cb||Vp|| i «- (i2) ||A 1 +f (w 2 )ll! 2(L2) + CxHA^t^ 2 )!! 2 ^ + Mo 

with < [3 < 1 and Co, C\, Mq depends only of the norm of initial data. 

As we have by energy inequalities ||Vp || £00(^2) < e < 1, we can conclude that: 

yp 2 \\ L 2 ( Hi+§) ^ M o 

with Mo depending only on the initial data. □ 

Corollary 1 Let the assumptions of theorem \ 3.1\ be satisfied with the exception of hy- 
pothesis on the viscosity coefficients which is replaced by: 
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it exists c > 0, sq > such that Vs > sq p(s) > c. 
it exists c\ > 0, m > 1 such that Vs > p(s) < 

i£ exists C2 > 0, m > 1 swc/i £/ia£ Vs > A(s) < - — . 



We get similar results as in theorem \3.1\ for all <p G C { 



WvfPW^fr+i) - M with - s < 2 ' 
where M depends only on the initial conditions data, on T, on cp, on s and on 
Proof: 

Indeed the only term which changes are the viscosity term. We have then in apply- 
ing the same operation as in the proof of theorem 13.11 the following terms to control: 

/ ipfi(p)D(u)A s (ipp 2 )dxdt and / / tpR ii j(\(p)d i Uj)A s ((pp 2 )dxdt. 
o Jr n Jo Jr n 

We conclude easily by the fact that y/ p(p)Vu and \J X(p)Vu belong to L 2 (L 2 ) and as 
p G L 00 ^ 1 ) we get that yW) and VKp) belongs to L°°(L P ) for all p > 2 by Sobolev 
embedding. We can then conclude by Holder inequality. □ 

3.2 Existence of global weak solutions for N = 2 away from vacuum 

We may now turn to our compactness result. First, we assume that a sequence (/?„, n n ) ng N 
of approximate weak solutions has been constructed by a mollifying process, which have 
suitable regularity to justify the formal estimates like the energy estimate and the pre- 
vious theorem. 

Moreover this sequence (p n ,u n ) n£ n has initial data ((po) n , (^o)n)) close to the energy 
space. In using the above energy inequalities, we assume that j-y({po)n), \^(po)n\ and 
(po)n\{uo)n\ 2 are bounded in L 1 ^^) so that (po) n is bounded in Ll(M. N ) . 

Then it follows from the energy inequality that: 

1. j 7 (p„),|V/9 n | 2 , p n \u n \ 2 are bounded in L°°(0, T, L 1 ^)), 

2. Du n is bounded in L 2 (R N x (0,T)), 

3. u n is bounded in L 2 (0, T, H 1 (Br)) for all R,T e (0,+oo). 

Extracting subsequences if necessary, we may assume that p n , u n converge weakly re- 
spectively in L 7 ((0,T) x B R ), L 2 {Q,T;H l {B R )) to p, u for all R,T G (0,+oo).We also 
extract subsequences for which ^/p~^u n , p n u n , p n u n (£> u n converge weakly as previously. 
Our goal now is to verify that the non-linear terms converge in the sense of the distribu- 
tion. The unique difficult term to treat is V p n ® Vp n . 

Moreover we assume that — is bounded in L 00 (L 00 ) and in using the previous theorem 
we get: 

Vy? G Co°(R iV ) (pp 2 n is bounded in L 2 (H 1+ ^) for all < s < 2. 
We can now show the following theorem: 
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Theorem 3.2 Let N = 2. We assume that there exists (3 > such that for all n G N: 

Pn(t,x) > (3 for a.a (t,x) G (0, +oo) x M 2 
Then there exists r] > such that if: 

\Np^ + \\VpSK\\\l^ + \\Mpo)\\l^<v 

then, up to a subsequence (p n ,u n ) converges strongly to a weak solution (p,u) (see def- 
inition \2.4\ ) of the system \1. Moreover we have Vp n (g> Vp n converges strongly in 
Lj oc (RxR N ). 

Proof of the theorem 13.21 : 

According to theorem I3T1 we have seen that for all <p G Cq°(R n ) <f(Pn — p 2 ) G L 2 {H 1+ i). 
We can now use some results of compactness to show that Vp 2 converge strongly in 
L 2 {L 2 oc ) to Vp 2 . We recall the following theorem from Aubin-Lions ( see Simon for 
general results [T5]). 

Lemma 1 Let X w Y > Z be Hilbert spaces such that the embedding from X in Y 
is compact. Let (/ n ) n gN a> sequence bounded in L q (0,T;Y), (with 1 < q < +oo) and 
(^ L ) n£ N bounded in L p (0,T;Z) (with 1 < p < +oo), then (/„) n eN ^ relatively compact 
in L q (0,T;Y). 

We need to localize because we have some result of compactness for the local Sobolev 

s_ 

space which is compactly embedded in Lf oc . Let (Xp)peN be a sequence of Cq°(R ) 
cut-off functions supported in the ball £?(0,p + l) of M. N and equal to 1 in a neighborhood 
of B(0,p). 

We have then in using mass equation: 

^y(Pn) + Vdiv(p 2 u n ) = -V(p 2 n dwu n ) 

We can then show that (^(Xp^Pn)))^^ is uniformly bounded for all p in L q {H a ) for 
a < in using energy inequalities and (XpV(p 2 )) ngN is uniformly bounded for all p in 

L 2 {Hi). Apply lemma ffl with the family (x P V(p 2 )) ngN and X = Xp H^ , Y = Xp L 2 , 
Z = XpH a , then use Cantor's diagonal process. This finally provides that: 

Vp > x P V(p2) ^ n _ +oc Xp V(p 2 ) in L 2 (L 2 ). (3.15) 

We now want to show some results of strong convergence for Vp n . 
We have then for all <f> G Cq°: 

\4>(V Pn -Vp)\ < —\ ( f > (Vp 2 n -Vp 2 )\ + \Vp 2 \\4>(—--)\=A n + B n . 

Pn Pn P 

We have A n converges to in L 2 (L 2 ) because y- is uniformly bounded in L 00 ^ 00 ) and in 
using (|3.15|) for all <f> G C^°, </>Vp 2 converges to Vp 2 in L 2 (L 2 ). In the same way we have 
— converges a.e to - (because we can show in using again the theorem of Aubin-Lions 
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that p n — p — > p — p in L 2 (L q ) up to an extraction with 1 < q < +00 then we can extract 
again a subsequence so that p n converges a.e to p). Moreover we have: 

I 1 1. 2 

' Pn P ~ P 

then by the theorem of dominated convergence we have B n tends to in L 2 (L 2 ). 
We can conclude that: 

V^GC °°, 4>Vp n ® V/> n 0V/> 8> Vp in L 1 (L 1 ). 

□ 



4 Existence of weak solution in the case N = 1 

We are now interested by the case N = 1. To start with, we focus on the gain of derivative 
for Vp 2 . 

4.1 Gain of derivative 

We can now write a theorem where we expose a gain of derivative on the density p in 
using the same type of inequalities as in the case N = 2. 

Theorem 4.3 Let (p, u) be a regular solution of the system with initial data in the 
energy space. Then we have: 

with < s < 5 and Mq depending only of the initial data. 

Remark 4 We observe the two important facts: 

1. We don't need any hypothesis on the size of the initial data. 

2. We don't need to localize because we know that p G Lf x . 

3. We don't need to assume that - G L°° , this is important. In fact if we assume 
that po admits some vacuum, we could show that Vp n (t,x) — > Vp a.e on the set 
A = {p > 0} because we have by compactness up to an extraction V/? 2 — ► Vp 2 a.e. 

Proof of theorem 14.31 : 

We use the same estimates as in the previous proof except for the delicate term: J K \d x p\ 2 A s p 2 . 
We have then d x p € L°°(L 2 ) and p - p G L°°(L 2 ) so p — p G L°°(F 1 ) and we have then 
by Sobolev embedding p G L 00 ^ 00 ). By composition theorem on the Sobolev space we 
get (p - p) 2 G L°°(H l ) and p 2 - p 2 = (p - p) 2 + 2p{p - p) G L 00 ^ 1 ). Finally we get for 
< s < 1, A s p 2 G L 00 ^ 1 - 3 ). Now for < s < \ by Sobolev embedding we obtain: 

AV G L 00 ^ 00 ). 
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So we can control the term j* J R \d x p\ 2 \A s p 2 as follows: 

f [ \d xP \ 2 \A s P 2 \ < lia^lll^). 

We treat the other terms similarly as in the previous proof. □ 
4.2 Results of compactness 

We can now prove our result of stability of solution in the case N = 1 in using the 
previous gain of derivative. Let (p n ,u n ) n eN a sequel of approximate weak solutions of 
system (II. ip . 

Theorem 4.4 Let (pq,Uq) initial data of the system lil.l]) in the energy space what it 
means that: 

' fQ\v*\* + h(Po) + \d x (Q\ 2 dx<M 



J 



with M > 0. 

Moreover we assume that p^ > c > 0. Then there exists a time T such that up to a 
subsequence, (p n ,u n ) converges strongly to a weak solution (p,u) on (0,T) x M in the 
sense of the distribution (see definition \2.J$ . 
Moreover d x p n converges strongly in L 2 (M x M.) to d x p. 

Proof of the theorem 14. 41 

We want now to control the vacuum of -. We recall that p- p G L°°(H l ) and §- t p 6 
L°°(L p ), then we have by Aubin-Lions theorem for all G C§°, <f>(p - p) G C([0, T],L°°). 
Moreover we set: 

fn{t) = / \d x p\ 2 lM.\B(0,n) 
JR 

then we have f n is a decreasing sequence and converges to 0. Then in using the theorem 
of Dini we get that on [0, T], f n converges uniformly to 0. Let e > 0, then it exists uq G N 
enough big such that for all n > uq: 

/ \9 X (P- P)| 2 lR\B(0,n) < e 
JR 

we proceed similarly to show that for a n' big enough we have for all n > n^: 

^|p-p| 2 l K \ B (0,n')^ £ 

then we have for G with a big enough support (p — p)(l — <j>) G L°°{H l ) with a 
norm inferior to e then it exists c such that by continuity with 0<e<c<p, \p — p\ < c. 
Finally we have ^ G L 00 (L 00 ) and we can conclude as in the previous theorem 14.31 show 
that: 

G C °° (Rp„ ^ n in L 2 (L 2 ). 

□ 
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Theorem 4.5 Let (pq,Uq) initial data of the system in the energy space. 

Then it exists e > such that if: 

\\d x pZh* + IIv^KHIls + \\jM)\\ L i < e 

then lip to (L SV,bs6C[U6TlCC (prio 

u n ) converges strongly to a weak solution (p,u) on R x R 

(see the definition \2.1$ . 

Moreover d x p n converges strongly to d x p in L 2 (R x R). 
Proof of theorem 14. 51 

We can show easily that \\p n — p||i,°o(#-i) < Ce and so we hav e — is uniformly bounded 
in L 00 (L 0C ). We can then conclude as in the proof of theorem 14.31 

□ 
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